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Abstract 

Measurements of the mean-square amplitude (9 2 ) of thermally-induced fluc- 
tuations of a thin layer of the nematic liquid crystal 152 subjected to a voltage 
V < V c are reported. They yield the limit of stability V c of the spatially uni- 
form conduction state to infinitesimal perturbations. Localized long-lived 
convecting structures known as worms form spontaneously from the fluctua- 
tions for V well below V c . This, as well as measurements of the lifetime of 
the conduction state for V < V c , suggests a thermally- activated mechanism 
of worm formation; but a theory for this non-potential system is lacking. 
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One of the uniquely nonlinear aspects of pattern-forming dissipative non-equilibrium 
systems is the occurrence of localized travelling-wave structures, or pulses, co-existing with 
the spatially uniform ground state. |I[ In one-dimensional cases |2|,[J it has been possible to 
understand such structures qualitatively in terms of a subcritical Ginzburg-Landau equation 

but for systems which are spatially extended in two dimensions || theoretical explana- 
tions are limited Recently, localized states were observed to co-exist with the ground 
state in a two-dimensional anisotropic system, namely in electro-convection (EC) of a thin 
layer of a nematic liquid crystal (NLC). HH These states, known as "worms", have a unique 
small width in one direction and a varying, usually much greater, length in the other. They 
are first observed as the control parameter of the system (the square V 2 of an applied volt- 
age) reaches a certain value; but initially they are very rare and become more abundant 
only as V 2 becomes larger. Their generation and decay as V was increased and decreased 
did not show any hysteresis. [§,[| Thus it was not possible to determine directly whether 
the creation of worms was associated with a supercritical or a subcritical bifurcation. A 
recently developed Swift-Hohenberg-like jit]] model equation [TTQ produced worm-like local- 



ized structures, but only for a strongly subcritical bifurcation. Worm-like solutions were 
found [JEJ also for coupled Ginzburg-Landau equations related to the equations of motion 
of the system (the weak-electrolyte model or WEM [0); but again the bifurcation to the 
non-uniform state was subcritical. 

In the present letter we report on a determination of the bifurcation point at which 
the spatially uniform state becomes unstable to infinitesimal perturbations by measuring 
the mean-square amplitude (8 2 ) of thermally-induced convection-roll fluctuations Jl4| about 



the non-convecting ground-state. This amplitude is known to diverge (except for nonlinear 
saturation) at the bifurcation point V 2 . A straight line through (6 2 )~ 2 vs. V 2 passes through 
zero at V 2 . We find that the worms grow spontaneously out of the ground state already for 
V 2 much less than V 2 , i.e. that they are associated with a subcritical bifurcation. We also 
measured the mean lifetime r of the conduction state as a function of the control parameter 
e = V 2 /V 2 — 1 and found that there is a wide range over which r decreases rapidly with 



increasing e. The data, as well as the spontaneous appearance of the worms at negative e 
and random spatial locations, suggest a thermally-activated mechanism for worm formation. 
However, the system under consideration is known to be non-potential, and thus it is not 
clear that the usual discussions of thermally-activated processes |15[ apply We believe that 
an explanation of our experimental data remains a challenge to the theory of non-equilibrium 
pattern-forming systems. 

NLC have an inherent orientational order, but no positional order. jl6| The average 



direction of their molecular alignment is called the director n. By confining a layer of NLC 
doped with ionic impurities between two properly treated glass plates, one obtains a cell 
with uniform planar alignment of the director parallel to the plates. If the NLC has a 
negative anisotropy of the dielectric constant and an AC voltage of amplitude V is applied 
between the plates, then a transition from the spatially uniform state to a convecting state 



with spatial variation occurs above a certain critical voltage V c . |17j The precise value of V c 
depends on the frequency / of the applied AC voltage and the conductivity a of the NLC. 
The concentration of ionic impurities determines a. 

The experiments reported here were done using the NLC J^-ethyl-2-fluoro-^ '-[2- (trans- 4- 
pentylcyclohexyl) ethylj '-biphenyl (152). Electroconvection in 152 leads to a great variety of 
spatio-temporal structures. [§,Pi|l8| Depending on a, one finds rolls, traveling waves, and 
spatio-temporal chaos at onset. The localized worm states mentioned above were observed 
when a was relatively small. 



The apparatus used for the experiment [18] consisted of a computer-controlled imaging 



system, temperature-control stage, and electronics for applying the AC voltage and mea- 
suring the conductivity of the cells. Only the component of the conductivity perpendicular 
to the director was determined. The thickness of the cells was d — 24 ± 2fim. The cells 
were uniform to about ±0.5/zm and the sample area was roughly 0.5cm x 0.5cm. Planar 
alignment was obtained by using a rubbed polyimide film which was spin-coated onto the 
transparent electrodes of indium-tin-oxide. The 152 was doped with molecular iodine. [T8 



The precise value of a was varied by changing the temperature. The conductivity decreased 



slowly over time typically at a rate of 1.5 x 10~ n (Ohm m day) -1 . This was compensated by 
slowly increasing the temperature, giving an almost constant conductivity over the duration 
of several months of the entire investigation. The measurements reported here were made 
for a = 3.7 x 10 -9 Ohm -1 m -1 at about 33°C and at 25Hz. 

The fluctuations and convection patterns were imaged using the shadowgraph technique 
19|,18| . An 8-Bit grey-scale frame-grabber was used to digitize the images. As the fluctua- 



tions are very weak, a series of 256 images ij(x, e) was taken at each value of e. The time 
between images was typically 5 s, and successive images were essentially uncorrelated. The 
images had to be taken before worms formed because the large worm amplitudes made it 
impossible to measure the fluctuation amplitudes. Close to onset worms formed very quickly. 
Thus for e > —0.08 it was impossible to accumulate enough images in their absence. For 
each image the signal 



7i(x,e) = 



7 i (x,e)-/ (x,e) // (x,e) (1) 



was calculated. Here Iq(x, e) is a background image obtained by averaging all 256 images. 
For each Ij(x, e) the structure factor (the square of the modulus of the Fourier transform) 
iSi(k, e) was calculated and the 256 Si(k, e) were averaged to get S(k, e). 

Figure [l] shows four examples of S(k, e) for different values of e. There are four peaks 



corresponding to two sets of rolls oriented obliquely to ft. p0| These two modes are known 



as zig and zag modes and are also seen in the extended-chaos state These examples 
show clearly that the peaks of S(k, e) get sharper as one gets closer to onset, as expected 
from theory ||21|| . They also get larger, but since the grey levels of each example in Fig. 1 
are scaled separately, this is not immediately apparent. 

In order to get the mean-square amplitude (A 2 ) of the signal fluctuations from the 
structure factor, one has to integrate over the peaks of 5(k, e) in k-space. The area under the 
peaks corresponds to (A 2 ). As the fluctuating pattern is very weak, there is a considerable 
background due to instrumental noise. We first calculated the azimuthal average S(k) from 
the data. Four examples of S(k) are shown in Fig. 0. We fitted a suitable function 



S(k) = {P{k) + B(k))/2nk (2) 

to this average. Here P(k) is the contribution from the peak and B(k) is the smooth 
background due to instrumental noise. A Lorentzian 

was used for the peak, and a polynomial in k was used for the background. The background 
included a term oc 1/k to account for the rapid increase of S(k) near k = noticeable in 
Fig. |2| (this increase is masked out in the examples in Figure p] in order to permit grey scales 
which make the fluctuations visible). The mean square amplitude (A 2 ) was then obtained 
from 

(a 2 }= r p(k) d k=^. (4) 



r 

For a comparison with theoretical results, (A 2 ) had to be related to the mean-square 
amplitudes (9 2 ) of the director-angle fluctuations. For our experimental setup (A 2 ) and {9 2 ) 
are related by ]T9" 

l + n)A n2 



(6 2 ) 



(A 2 ) . (5) 



And 

Here h — 1 — {n e /no) 2 with n e and no the index of refraction parallel and perpendicular to 
the director, d is the cell thickness, and A is the wavelength of the pattern. 

Figure |3] shows (6 2 )~ 2 vs. V 2 . As expected, the data fall on a straight line. A fit of a 
straight line to the data gives V 2 = 663 Volt 2 or V c = 25.7 Volt. Using this value to compute 
e, we get the results shown in the top half of Fig. |4j. A straight line through the data then 
corresponds to (9 2 ) = 6 2 e~ 1 ' 2 , with #0 = 3.07mrad. A simple theoretical estimate for Q Q is 
given by 

<"»> - k £ d < 6 > 

where fc& is Boltzmann's constant and k e i = 20.7 x 10™ 12 iV an average orientational elasticity 



of 152 [rS|. This gives 6$ = 2.9 mrad, which is in good agreement with our results. This 



agreement confirms the assertion that the fluctuations are of thermal origin. 



From P(k) it is possible also to extract an average of the two-point correlation length. 
It is given in terms of the half width r by £ = 1/T. Theory predicts |I| that the dependence 
of the correlation length on e should be £ = £o( —e ) The bottom half of Fig. ^ shows 
l/£ 2 vs e. The data points fall on a straight line through e = 0, consistent with the expected 
behavior. The amplitude £o is found to be 0.31<i. 

We found that worms form spontaneously from the fluctuations of the conduction state 
already well below V c . Thus, the bifurcation associated with them must be subcritical. To 
get more insight into the nature of worm formation, we measured the mean lifetime r of 
the conduction state as a function of e. A voltage was applied to the cell and then the time 
until a worm first appeared was measured. This was repeated many times and at various 
voltages to get a statistical average for r(e) [p^] . Figure || shows the frequency 1/r, in 
p3fl These results are reminiscent of thermally-activated nucleation in potential systems. 
Thus we show in Fig. [| a simplified version of Kramers' formula [|15j for the diffusion of a 
particle in a two- well potential, namely 1/r = f Q exp(—AF). Here AF is the value of the 
Landau potential F = — \eA 2 + \gA^ + \kA & at the unstable fixed point which is located 
at A 2 = (—g — \/ g 2 + Ake)/2k (the effective noise strength is absorbed in the scale of F). 
Guided by the data, we chose the saddlenode of F to be located at e = —0.5, and used 
k = 0.00034 and fo = 0.0088 s _1 . This model reproduces the trend of the experimental 
results remarkably well. The deviations near e = where AF becomes small are expected. 
We do not feel that a more sophisticated analysis in terms of a potential model would be 
justified. Instead, a theory for non-potential systems should be developed if possible. 

We have shown that measurements of the critical behavior of thermally-induced fluctu- 
ations can yield the critical voltage V c of the conduction state. We found that the director- 
angle fluctuations (8 2 ) agree quantitatively with the result expected for thermally-driven 
fluctuations. Worms appeared already well below V c , indicating that the bifurcation from 
conduction to worms at low conductivities is strongly subcritical. At higher conductivities a 
supercritical bifurcation had been found to an extended-chaos state || . An interesting ques- 
tion is how the bifurcation evolves from sub- to super-critical as a increases. We determined 
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also an angular average £ (e) of the two-point correlation length. It agrees with the expected 
relation £ = ^e' 1 ^ 2 and yields £ = 0.31. We measured the lifetime of the conduction state, 
and found it to be finite for e > —0.6 and rapidly decreasing as e approached zero from 
below. The results are reminiscent of thermally activated nucleation. Surprisingly they can 
be fit by a potential model even for this non-potential system where there is no firmly based 
theory of which we are aware. 

This work was supported by the National Science Foundation through grant DMR94- 
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FIGURES 




FIG. 1. Time-averaged structure factors S(k) of the fluctuations for (a) e = —0.46, (b) 
e = —0.21, (c) e = —0.11, and (d) e = —0.08. The four peaks correspond to the zig and zag 
modes which are seen also in the extended-chaos state for large conductivity and e > 0. The 
director is in the horizontal direction. 
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FIG. 2. The structure factor S(k) for e = -0.08 (solid squares), -0.11 (open squares), -0.21 (solid 
circles), and -0.46 (open circles). 



11 




FIG. 3. The inverse square of the mean-square director-angle fluctuations < 9 2 >~ 2 as a 
function of the applied voltage- amplitude V 2 . The intercept of the straight-line fit to the data 
with the horizontal axis gives V 2 . 
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FIG. 4. Top: The inverse square of the mean-square director-angle fluctuation < 2 > 2 as a 
function of e = V 2 /V 2 - 1. Bottom: The inverse square l/£ 2 of the average two-point correlation 
length derived from the width of S(k) in Fig. [||. 
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FIG. 5. The inverse lifetime 1/r of the spatially uniform state vs. e. The solid line is based on 
a potential thermal-activation model. 
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